Lower Bound on the Magnetic Field Strength of 
a Magnetar from Analysis of SGR Giant Flares 

- ; A. A. Gvozdevt I. S. Ognev J E. V. Osokina 

^ I Yaroslavl State University, Sovietskaya 14, 

(N ! 150000 Yaroslavl, Russia. 

>>: 

O . Abstract 
(N ■ 

! Based on the magnetar model, we have studied in detail the processes of neutrino 

' cooling of an electron-positron plasma generating an SGR giant flare and the in- 

jnr^ . fluence of the magnetar magnetic field on these processes. Electron-positron pair 

^ I annihilation and synchrotron neutrino emission are shown to make a dominant con- 

Q^' tribution to the neutrino emissivity of such a plasma. We have calculated the 

Q ■ neutrino energy losses from a plasma-filled region at the long tail stage of the SGR 

'. 0526-66, SGR 1806-20, and SGR 1900+14 giant flares. This plasma can emit the 

I energy observed in an SGR giant flare only in the presence of a strong magnetic 
field suppressing its neutrino energy losses. We have obtained a lower bound on 

^ . the magnetic field strength and showed this value to be higher than the upper limit 

^ I following from an estimate of the magnetic dipole losses for the magnetars being 

■ analyzed in a wide range of magnetar model parameters. Thus, it is problematic to 
. explain the observed energy release at the long tail stage of an SGR giant flare in 

Ti^lj- I terms of the magnetar model. 

in : 
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> ■ INTRODUCTION 

^ ■ Soft gamma-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs) constitute a 

. special class of neutron stars with anomalously large spin periods P ~ (5 — 8)s and 

spindown rates P ~ (10"^'' — 10~^^)ss~^. About 20 such objects have been discovered to 
date in our and nearest galaxies (Mereghetti 2008). Their identification with supernova 
remnants shows them to be young isolated neutron stars with ages tns ~ (10^ — 10^) yr 
without accretion disks (Bisnovatyi-Kogan 2006). Note the most characteristic properties 
of these objects. First, in quiescence they have an anomalously high effective temperature 
for isolated stars, emitting soft X rays from the surface with luminosities L^r^ ~ (10^^ — 
10^^) (Mereghetti 2008). Second, AXPs exhibit large glitches (Dib et al. 2008), whileSGRs 
exhibit gamma-ray bursts (Strohmayer and Watts 2006; Watts and Strohmayer 2007). 
They are interpreted as a manifestation of a seismic activity in these objects similar to 
the seismicity of the Earth, the Sun, and the young Vela pulsar (Gogus et al. 2000). Third, 
gamma-ray bursts were also detected from AXPs, although they are not so powerful as 
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those from SGRs (Woods et al. 2005). This suggests that SGRs and AXPs most hkely 
belong to the same class of neutron stars. The SGR flare activity manifests itself in the 
emission of numerous (up to 100 episodes per day) short bursts in the energy range from 
hard X rays to soft gamma rays with a typical duration of Tp ~ 0.1 and energy Ep < 10^^ 
erg. In several cases, a series of short bursts was followed by a giant flare exceeding 
in energetics the short ones by several orders of magnitude. Two easily distinguishable 
stages were observed in the three most powerful (in energy release) giant flares from SGR 
0526-66 (March 5, 1979), SGR 1900+14 (August 27, 1998), and SGR 1806-20 (December 
27, 2004): a short, with a duration ths ~ (0.25 - 0.5) and energy Ehs ^ (10^^ - 10*^) 
erg, hard spike (HS) followed by a pulsating long (t^t ~ (200 — 400) )s tail (LT) with 
energy Ei^ ~ (1 — 4) x 10*^^ erg, at which a modulation of the emission intensity by the 
neutron star spin period was observed. Here, we investigate the energy losses at the LT 
stage of an SGR giant flare and do not consider the HS stage. Below, we present data for 
the three most energetic SGR giant flares at this stage: 

SGR 0526-66 {D ^ 55 kpc) tlt ~ 200 s, Elt ~ 3.6 x 10^^ erg (Mazets et al. 1979); 

SGR 1900+14 {D ^ 15 kpc) tlt ~ 400 s, Elt ~ 1-2 x 10^^ erg (Ibrahim et al. 2001); 
SGR 1806-20 {D ^ 15 kpc) tlt ~ 380 s, Elt ~ 1.3 x 10^^ erg (Mereghetti et al. 2005; 
Frederiks et al. 2007). As can be seen, these parameters almost coincide, suggesting a 
unifled flare formation mechanism at the LT stage. 

Note that the energy release in SGR giant flares is smaller than that in supernova 
explosions and cosmological gamma-ray bursts only. The model of a magnetar, a neutron 
star with an anomalously strong magnetic fleld Bm ~ 10^^ G (Duncan and Thompson 
1992; Thompson and Duncan 1993), was proposed to explain such a huge gamma-ray 
burst energy for an isolated neutron star. In this model, the energy of the magnetic 
field liberated when its configuration changed rapidly is assumed to be the gamma-ray 
energy source. Subsequently, the magnetar model was used to describe the emission 
from AXPs and SGRs in quiescence (Thompson and Duncan 1996) and during an SGR 
giant flare (Thompson and Duncan 1995). In the magnetar model of an SGR giant flare 
(Thompson and Duncan 1995), the long-term {tns ~ lO^yr) evolution of the magnetar 
core with poloidal and toroidal magnetic fields of strength Bm ~ 10^^ G is assumed 
to be ended with a starquake leading to a largescale plastic deformation of its crust. 
The electric currents emerging during the deformation produce a perturbative magnetic 
fleld with fleld lines closed on the crust. The region with closed fleld lines is rapidly (in 
hundredths of a second) flUed with an electron- positron plasma trapped by this fleld (the 
so-called fireball). A fairly hot plasma, with a temperature T ~ 10^*^ K, is generated. The 
X-ray photon flux observed at the LT stage is assumed to be emitted from a thin near- 
surface layer of the flreball. Thompson and Duncan (2001) modeled this emission and 
compared it with the observed light curve of the SGR 1900+14 giant flare. As a result, 
the distributions of plasma parameters (temperature and magnetic fleld) that agreed best 
with the observational data were obtained. It is important to note that the neutrino 
emission was disregarded by these authors, because it was assumed to be significantly 
suppressed by a strong magnetic fleld inside the flreball. 

Here, we study in detail the plasma neutrino emission processes at the LT stage based 
on the magnetar model of an SGR giant flare. We show that the plasma energy losses 
through neutrino emission are significant even in the case of a strong magnetic field with 
a strength B > 10^^ G. For the most energetic fiares from SGR 0526-66, SGR 1806-20, 
and SGR 1900+14, we model the flreball neutrino cooling. The dependences of neutrino 
cooling on parameters of the temperature and magnetic field distributions, fireball size, 
and its total energy are analyzed. 
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Below, except for the specially stipulated cases, we use a system of units in which 
c = h = k = 1. 



NEUTRINO COOLING OF A RELATIVISTIC 
NONDEGENERATE ELECTRON-POSITRON PLASMA 

Main Neutrino Processes 

The following reactions are the most significant neutrino emission processes of a rela- 
tivistic nondegenerate electron-positron plasma: the electron-positron pair annihilation 
into a pair of neutrinos with an arbitrary flavor 

e" + e+ -)■ z/j + Ui, (1) 

the plasmon decay into a neutrino pair 

7 z/j + z/i, (2) 

the neutrino production due to the fusion of two photons 

7 + 7 z/i + z/i, (3) 

the photoneutrino emission process 

+ 7 ^ + z/i + Di, (4) 

the neutrino synchrotron emission by electrons 

A + z/j + z>j, (5) 

which is kinematically allowed only in an external magnetic field. Here and below, the 
subscript z = e, /i, r specifies the neutrino flavor. 

The Case of a Weak Magnetic Field 

We will begin our analysis of the neutrino emission processes with the case of an 
electron-positron plasma in the absence of a magnetic field. The neutrino emissivity of 
an ultrarelativistic nondegenerate plasma in the annihilation process (1) is well known 
(Kaminker et al. 1992) and can be written as 

Qio) = im^GlT^. (6) 

iZ 71 

Here, T -is the plasma temperature, Gp -is the Fermi constant, = J2i (^^^ + '^aj — 
1.675, where c„. and Ca. are the vector and axial constants of the leptonic electroweak 
current, and is the Riemann zeta function. 

The emissivity in the plasmon decay into a neutrino pair (j2]) has also been well studied 
and can be found, for example, in Yakovlev et al. (2001) and Kantor and Gusakov (2007). 
For a nondegenerate plasma, it can be represented as 

QP = ^a^GlT^ (7) 
o24 IT 
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where a = 1/137 -is the fine-structure constant and = J^i'^li = 0.9248. It is easy to 
see that the emissivity in this process is strongly suppressed compared to the electron- 
positron pair annihilation. 

The processes ([3]) and (jl]) in a plasma were extensively studied previously (Beaudet et 
al. 1967; Itoh et al. 1996) and the neutrino emissivities in these processes were shown to 
be also negligibly small compared to (6). Thus, the annihilation process makes a major 
contribution to the neutrino emissivity of a relativistic nondegenerate plasma. 

In the case of a relatively weak magnetic field, ^ eB the emissivities in the 
processes ([I]), ([2]), ([3]) and (jl]) change insignificantly, but the presence of a magnetic field 
makes the new synchrotron neutrino pair production process ([5]) kinematically open. The 
emissivity in this process in the limit of a weak magnetic field is given by the expression 
(Kaminker and Yakovlev 1993): 



^(0) _ 10C(5) ^2 ^2 / Tj\2 rp5 



In (11) +4.66 



from which it follows that Q's'' IQ^a ~ {eB/T'^Y i-^-, the emissivity in this process in the 
limit under consideration is also suppressed. Thus, in a relatively weak magnetic field, the 

pair annihilation reaction is the main plasma neutrino cooling process and the total 
neutrino emissivity is defined by Eq. 

The above analysis allows the characteristic plasma neutrino cooling time to a 
temperature T to be calculated. It can be found from the equation 

"llTT^ 



d 



60 



-T 



Qf. (9) 



where the total energy density of an ultrarelativistic nondegenerate plasma and photons 
appears on the left-hand side. The solution of this equation, which essentially coincides 
with the estimate from Thompson and Duncan (1995) 

m 44 TT^ 1 _ /lMeV\^ 



shows that the neutrino cooling time for a fairly weak magnetic field is an order of mag- 
nitude shorter than the characteristic duration of a giant fiare tlt — (200 — 400) s. Thus, 
the bulk of the hot-plasma energy is expended in cooling by neutrino emission and such a 
plasma cannot be the source of an SGR giant fiare. Consequently, a mechanism suppress- 
ing the neutrino cooling processes is required. As will be shown below, a strong magnetic 
field that is capable of significantly reducing the hot-plasma energy losses through neu- 
trino emission can act as such a mechanism. 



The Case of a Strong Magnetic Field 

In the asymptotic limit of a strong magnetic field, when eB ^ > m? , the emissivity 
in the electron-positron pair annihilation process (1) is well known and can be represented 
for a nondegenerate plasma as (Kaminker et al. 1992) 

= ^§^G|m^e5T^ (11) 

where m is the electron mass. 
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The plasmon decay reaction is also modified significantly by a strong magnetic field, 
because not only the amplitude of the process but also the dispersion law of plasmon modes 
change. Calculations show (Kuznetsov et al. 1998) that the emissivity in the strong-field 
limit is 

Q^P = ^^aGl{eBfT\ (12) 

As we see from this expression, the emissivity in this process can become equal to (11) 
only at a sufficiently large magnetic field strength i? ~ 8 x 10^^ G. 

The fusion of two photons ([3]) in the case of a strong magnetic field was investigated 
by Rumyantsev and Chistyakov (2008). Since the analytical expression for the neutrino 
emissivity in this reaction is fairly cumbersome, it is given below in an approximate form 
that is valid only in the limit eB ^ ^ m?: 

Q(^)^2.7xl0^«^f-y. (13) 
' ' cm-* s \m J 

Note that the emissivity in this limit is virtually independent of the magnetic field strength 
and is lower than that in the annihilation process approximately by three orders of mag- 
nitude. 

The emissivity in the photoneutrino production process (jl]) in the limit of a strong 
magnetic field can be represented as 

2/72 

Q^^^ c:, —talGlm^ eB In (T/m) , (14) 
where the numerical factor I has the following integral representation: 

au av X 

, XU _j_ g XV 



Jo Jo 



X |(1 - t;) [(1 - vf - (1 - uf] e^" + {l + v) [{1 + vf - (1 - uf] 



e-^n -0.09. 



Comparison of the derived emissivity with ffTTj) shows a negligible contribution from this 
reaction to the neutrino energy losses of a strongly magnetized nondegenerate plasma. 

The neutrino emissivity in the synchrotron emission process in the strong-field 
limit was obtained by Kaminker and Yakovlev (1993): 

Qs'' = ^r^ClGl{eBfy^T'/^e-^/^. (15) 

As we see from this expression, the emissivity in this limit is exponentially suppressed by 
the smallness of the number density of electrons and positrons at all Landau levels, except 
for the ground one. On these grounds, it is generally concluded that the synchrotron 
emission process cannot play a significant role in the neutrino cooling of a plasma with a 
strong magnetic field. However, as will be shown below, this conclusion is unjustified. 

Assuming that the electron-positron pair annihilation ([1]), is the dominant plasma 
neutrino cooling process in the presence of a strong magnetic field, we can find the charac- 
teristic neutrino cooling time r!) to a temperature T. It can be found from the equation 
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where the energy density of an electron-positron plasma whose particles are at the ground 
Landau level appears on the left-hand side. Solving this equation gives 



Note that the neutrino cooling time depends only on temperature, because both the 
plasma energy density and the energy losses through neutrino emission in this limit grow 
proportionally to the magnetic field strength. It follows from the above estimate that 
the characteristic neutrino cooling time for a plasma with a strong magnetic field must 
exceed the giant flare duration ~ (200 — 400) s by several times. On these grounds, 
the authors of the magnetar model concluded that the plasma neutrino emission could be 
neglected (Thompson and Duncan 1995). 

Although the neutrino cooling time in this limit does not depend on magnetic field 
strength, the neutrino emissivity ( ITTl) increases linearly with growing field, reaching the 
field-free value ([6]) at i? ~ 10^^ G {T/mY. Thus, the neutrino emissivity in the annihilation 
process must be suppressed significantly in a magnetic field whose strength satisfies the 
following inequality: 

4.4 X lO^^G < < lO^^^G or ^2 < x < 22 t, (18) 



where t = T/m, x = y/2eB/T. The lower bound follows from the condition eB ^ T^. 
Note that the commonly assumed magnetic field strength for magnetars Bm ~ 10^^ G 
(Duncan and Thompson 1992; Thompson and Duncan 1993) falls within this range at a 
plasma temperature T > m typical of a giant flare. However, the estimate ( !T7|) for the 
neutrino cooling time is valid only if the neutrino emissivity is defined by the asymptotic 
expression ( ITTj) . A detailed analysis shows that this expression is applicable only in 
magnetic fields with strengths 

B> 7 xlO^^ Gt^, (19) 

while the emissivity in this process at smaller strengths can exceed considerably the 
asymptotic one. 

In addition, the neutrino synchrotron emission process ([5]) is important for plasma 
cooling even in a strong magnetic field. Indeed, the ratio of the emissivity in this process 
to the emissivity in the annihilation reaction is given by the expression 

_ ^/2^ fi-l)t' x''/' e-^ (20) 



Q(f) 37r2C(3) V 4e 

which has a maximum at Xmax = 13/2 equal to Q^P /Q^a'' — 3-5 ■ Consequently, at 
a temperature T > m typical of a giant flare, the neutrino synchrotron emission pro- 
cess contributes significantly to plasma neutrino cooling in the range of magnetic field 
strengths f lTS]) . where the annihilation process is essentially suppressed. 

Thus, our analysis shows that to describe the neutrino cooling of a plasma emitting 
an SGR giant flare, it is insufficient to consider only the electron-positron pair annihi- 
lation ([1]), because the synchrotron neutrino pair production ([5]) is a no less important 
reaction and both these processes make a comparable contribution to the plasma energy 
losses. In addition, the actual neutrino emissivities in the processes ([T]) and (E]) exceed 
considerably the asymptotic emissivity fllip in the strong-field limit even in magnetic 
fields with strengths B > 10^^. 
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LOWER BOUND ON THE MAGNETIC FIELD STRENGTH 
OF A MAGNETAR FROM THE NEUTRINO COOLING RATE 

In this section, the resuhs obtained above are used to model the plasma neutrino 
cooling for an SGR giant flare. First of all, it should be emphasized that, although such a 
plasma is fairly hot, it is essentially transparent to neutrinos. Consequently, the neutrinos 
escape freely from the entire plasma-occupied volume and the energy lost through neutrino 
emission does not depend on the geometry of the emitting region. To simplify our analysis, 
we model the fireball as part of a sphere with radius Rq whose center is on the magnetar 
surface (Fig. [1]). The distributions of plasma parameters are assumed to be spherically 
symmetric, just as in the paper by Thompson and Duncan (2001). In our numerical 
calculations, we used the following temperature and magnetic field strength distributions 
inside the fireball: 

t{z)=to{l + z)\ (21) 

6(^) = 6o(l + ^)^ (22) 

where z = r / Rq is the distance from the fireball center in units of its radius, t = T/m and 
b = eB/m? are the dimensionless temperature and magnetic field strength, respectively. 
The parameters of the distributions tQ,bo,P and 7 allow the neutrino emission to be 
described completely. However, in general, we cannot use the results from Thompson 
and Duncan (2001), where these parameters were obtained by comparing the modeled 
X-ray emission with the observations of a giant flare from SGR 1900+14, because the 
neutrino emission was completely excluded in this analysis. Indeed, the plasma at such 
parameters has only the energy that was observed in this flare as X-ray emission. In 
contrast, in the case where the energy losses of the medium through neutrino emission are 
important, the plasma energy must be considerably higher. Thus, the parameters of the 
distributions (12T|) and (122|) must be found by simultaneously modeling the X-ray emission 
from the fireball surface and the neutrino cooling of the medium from its volume and by 
comparing the modeling results with the observed X-ray light curves of SGR giant flares. 

Here, we consider a simpler model. The electron-positron plasma is assumed to have 
an energy Etot that is expended in cooling by neutrino emission and an SGR X-ray flare 
with an energy Elt- In this case, the energy balance equations can be written as 

1 

Etot = vElt = 2nRl j {u,±{z) + f/^(z)) z^ dz, (23) 



1 

E, = iri- 1) Elt = 2tiRI tlt j Qu{z) z" dz, (24) 



where 77 = Etot/EiT the plasma is assumed to be composed of electrons, positrons, and 
photons, Ue±{z) and U^{z) are the local energy densities of these particles, Qy{z) is the 
neutrino emissivity, and the right-hand side of the second equation gives the total energy 
loss Ey through neutrino cooling in the giant-flare time tlt- 

Analytical Neutrino Emission Model in the Asymptotic Limit of a Strong 
Magnetic Field 

As was shown above, the asymptotic expressions for the neutrino emissivities derived 
in the limit of a strong magnetic field cannot be used to describe the neutrino emission 
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of an SGR giant flare. However, we will consider this limiting case, because it admits an 
analytical solution and provides an insight into the main features of neutrino cooling in 
the magnetar model. In this case, the neutrino energy losses are determined only by the 
electron-positron pair annihilation ([1]): 

Qu = Q^a\ (25) 

whose emissivity is given by Eq. (11). The plasma energy density in this limit is deter- 
mined by electrons and positrons at the ground Landau level: 

4 

U,± + U^c^Ue±-—bf. (26) 
The energy balance equations (|23|) and fl2^ in this case can be represented as 



vElt = 1 J{f3 + 27) m^Rl bo tl, (27) 
if] - 1) Elt = Af3 + 57) Gl TLT Rl bo tl (28) 



where the function J{6) is defined by the integral 



and the previously introduced temperature (1211) and magnetic field strength 0221) distri- 
butions are used. The solution of this system of equations can be represented as 

- ymCl J(/3 + 57) G| T,r ^) ' 

^ 3(C(3)Ci)^/^J^/3(/j-^57) gfE..r;/? 
"ol'/J 21/3 j5/3(/3 + 27) m2/3i?3 (ry- 1)2/3- ^''"J 

Analysis of this solution shows that to(^) and (3 depend very weakly on the exponents bo{rj) 
and 7, respectively. In addition, in the magnetar model (Thompson and Duncan 2001), the 
change in temperature inside the fireball is assumed to be fairly small from the isothermal 
case (7 = 0) to the case of 7 = —1. Of particular interest is 7 = —1/2 corresponding to 
an arbitrary exponent f3. Since a sharp decrease in magnetic field strength in the fireball 
would lead to excessively intense cooling of its outer layers through neutrino emission, 
we will assume that the field decreases no faster than the dipolar law /3 = —3. Under 
these assumptions, the following approximate expressions can be used to estimate the 
temperature To ~ 0.51 MeV ■ to and magnetic field strength Bo — AA x 10^^ G - bo aX the 
fireball center: 

1 / — 1 \ 

to(r7) ^4.4(1- 0.67) ^ , (31) 

''"100 V ^ / 

bo{^) ^ 2.1 (1 - 0.5/3 + 0.3/3^) ^ (32) 

where xioo = Ti^/lOO s, £"44 = E'^t/IO^^ erg and 10 = -Ro/10 km. As we see from these 
expressions, the solution toirf) grows rapidly with plasma energy Etot = V ^lt in the 
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case of relatively small neutrino energy losses, when r/ ^ 1. Thus, small neutrino energy 
losses are possible only in a fairly cold plasma. In the region where neutrino cooling 
dominates, when ^ 1, the solution toif]) reaches a constant. The existence of this 
limiting temperature stems from the fact that the change in fireball temperature through 
the emission of neutrinos and photons in the flare time tlt is neglected in the model 
under consideration. The solution &o(^) contains a divergence at 77 = 1 which physically 
corresponds to the impossibility of completely removing the energy losses through neutrino 
emission. It passes through a minimum at rj = 5/3 and reaches the asymptotics rj as 
boir]) ~ rj increases further in the range of energies where the neutrino energy losses 



dominate. Figure [2] presents the corresponding solutions at 7 = —1/2 and /3 = —3 for 
the following giant flares: 

SGR 0526 - 66 : E44 = 3.6, rioo = 2.0, (33) 

SGR 1806 - 20 : E44 = 1.3, noo = 3.8. (34) 



We do not discuss the flare from SGR 1900+14, because it has characteristics at the LT 
state similar to those of SGR 1806-20. 

Let us separately consider another important peculiarity of the solution obtained. If we 
fix the total plasma energy Etot = V^lt and assume that the temperature to > to{ri) then, 
in view of Eq. fl27j) . the magnetic field strength Bq must be lower than the corresponding 
solution bo{ri). It is easy to see that, in this case, the energy losses by such a medium 
Eu through neutrino emission that are defined by the right-hand side of Eq. (125]) will 
increase, i.e., they will become greater than {Efot — Elt)- In contrast, for to < ^o(^) the 
value of 60 must be higher than 6o('?) and the neutrino energy losses will decrease. Thus, 
the solution (129|) defines the maximum temperature, while the solution (130|) defines the 
minimum magnetic field strength at which the neutrino energy losses leave an energy in 
the plasma no less than E^t observed in the photon emission. In addition, since to(^) 
has a global maximum and 60 (^7) has a global minimum, the upper and lower bounds, 
respectively, on the admissible temperature and the magnetic field strength at the fireball 
center correspond to them: 

y(max) ^ MeV (1 - 0.67) (35) 

'100 

Sj™") ~ 2.8 X 10^^ G (1 - 0.5 /3 + 0.3 Z^^) =^^iji^, (36) 

They are needed for the energy emitted by the plasma in photons to be no less than Elt- 
As was noted above, the existence of a maximum temperature of the medium stems from 
the fact that its change through the emission of photons and neutrinos in the flare time 
tlt is disregarded in the model under consideration. It follows from the condition ( 136|) 
that, depending on the magnetic field strength, a plasma with an arbitrary energy Etot 
can emit an energy in photons no greater than 

^(^a.)^^Q44g^^:Bi^^ (37) 
'^100 

where -B15 = Bq/IQ^^G and /3 = —3. As we see from this expression, the fireball radius Rq 
is the most important parameter defining the admissible energy release in photons, while 
the reasonable energy release itself is close to the typical energy of SGR giant flares. 
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It should be noted that a modulation of the emission intensity by the neutron star 
spin period (Mazets et al. 1979; Ibrahim et al. 2001; Mereghetti et al. 2005; Frederiks et 
al. 2007) is clearly traceable for all of the known SGR giant flares at the LT stage. The 
presence of such pulsations leads to the conclusion that the size Rq of the plasma-occupied 
region must be close to the neutron star radius Rns- Under this assumption and for a 
standard value of R^s — 10 km and a typical energy E^t ^ 10^^ erg of an SGR giant 
flare at the LT stage, the magnetic field strength must be Bq > 10^^ G. It is interesting 
to compare this value with the upper bound on the magnetic field strength that follows 
from an estimate of the magnetic dipole losses for a magnetar: 

1 M"*^^^ / ■ \ 1/2 
Sa/d=^2.1x101^G'--^ PioP-io . (38) 

Here, Mi,4 = MmsI^-^Mq, Rw = RNs/^^km, Pio = P/lOs, P_io = P/IO^^^, where 
Mj^Si P-i and P are the neutron star mass, period, and spindown rate, respectively, 6 is 
the angle between the angular velocity and magnetic moment vectors. It is easy to see 
that these bounds almost coincide at Rq ~ Rns- Thus, the above estimates leave open 
only a narrow range of magnetar magnetic field strengths even in the simplified model 
that grossly underestimates the energy losses of the medium through neutrino emission 
and, consequently, the minimum magnetic field strength b'^^^\ 



Modeling the Neutrino Cooling of SGR 0526-66 and SGR 1806-20 

In the previous section, we studied an analytical model for the neutrino cooling of a 
plasma producing an SGR giant flare that is based on the asymptotic expressions for the 
plasma energy density and neutrino emissivity derived in the limit of very strong magnetic 
fields. In this section, we will use the expressions valid for an arbitrary magnetic field 
strength for a more realistic modeling of the neutrino cooling process. 

Since the photon mean free path in the plasma under consideration is small, the 
electrons, positrons, and photons are in local thermodynamic equilibrium and, hence, the 
equilibrium distribution functions can be used to describe them. In this case, the plasma 
energy density will be defined by the expressions 



4 

m 



TT m 



15 



oo oo 

J eV^A + 1 



X 



+ 26n+ 1)1/2 



QVx'^+2bn+l/t _|_ ]_ 



dx 



(39) 
(40) 



As was shown above, the electron-positron pair annihilation ([T]) and the neutrino syn- 
chrotron emission (|5]) are the main neutrino cooling processes for the plasma under 
consideration: 

Qu = Qa + Qs- (41) 



Here, for the neutrino emissivity of the annihilation process Qa, we used the interpolation 
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formula (Kaminker et al. 1992) 



15 



t]t^ + Cl P{t) 



F{t,b) 



(42) 



+ 



{Cl - Cl) - 



62 Sit) 

mrT6) r^n-t 



Pit) =t^{l + 3.581 1 + 39.64 1^ + 24.43 + 36.49 + 18.75 1^) 
S{t) = t (l + 1.058 1 + 0.6701 1^ + 0.9143 1^ + 0.472 1^) , 



F{t,h) 



1 



■, i?i = 1 + Q 



exp 



26 



3t 



where Cl ~ 1.675, ~ 0.175 the constants ci ~ 3.106 x 10"^ C2 ~ 1.491 x 10-^ 
C3 ~ 4.839 X 10^^. For the emissivity of the synchrotron neutrino pair production reaction, 
we used the interpolation expression from Kaminker and Yakovlev (1993) 



6^ 



;i + ai6)= 



02 t 



G|m9 
120 7r3 

2 



N,± {ClF+ - CIfJ) 



(43) 



n=l 



1 + 



dx 



a^tb 



:i + ai6)V2; ■ (1 

^ 6^(1 + get) 
(l + ai6)6 



1 + 



ai6)3/2 

a^tb 



1 + In 1 + 
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where is the number density of electrons and positrons at all Landau levels, except 
for the ground one. Here, the constants Oi ~ 0.955, 02 — 9.439, 03 ^ 23.31, ~ 0.26, 
ag ~ 0.168, ae ~ 0.971. 

We modeled neutrino cooling for giant flares from SGR 0526-66 and SGR 1806-20 
with the characteristics ( l33i) and ( IMI) . As was noted above, the SGR 1900+14 flare at 
the LT stage is similar in its properties to SGR 1806-20 and below we do not separate 
them. Our numerical calculations confirm all of the qualitative conclusions reached in the 
previous section for the analytical model of neutrino cooling, but the numerical values 
of the plasma temperature and magnetic field strength change significantly. Below, we 
present the results for the case of /3 = — 3 which corresponds to a dipole magnetic field 
configuration in the fireball. The law of change in temperature was chosen with an 
exponent characteristic of the magnetar model, 7 = —1/2 which was compared with 
the isothermal case of 7 = 0. The fireball radius was chosen to be Rq = 10 km, which 
corresponds to the standard radius of a neutron star. The results of our numerical solution 
of the system of equations fl23l i^^ for the SGR giant flares under consideration are 
presented in Figs. [3] and IH As we see from the plots, the general trend in the behavior 
of the temperature and magnetic field strength as a function of the total plasma energy 
Etot = V^LT is the same as that in the analytical model. However, the plasma being 
analyzed turns out to be colder and the magnetic field strength in it must be higher by 
several times. Note that the dependence of the solutions on parameters 7 and (3 also 
remains similar to the analytical model (see Eqs. fl3ip and (1321) ). Thus, the calculated 
magnetic field strength in the fireball is virtually independent of the law of change in 



11 



temperature, as is demonstrated in Fig. |H We emphasize once again that the derived 
field strength at fixed total plasma energy is minimally possible for the emission of the 
observed giant-fiare energy Elt- Note also that the global minimum of the magnetic field 
strength for the SGR giant fiares in question takes place at ~ 1.3 instead of 77 = 5/3 in 
the analytical model and, hence, the energy losses through neutrino emission are reduced 
approximately by half. However, the losses remain fairly large and their reduction requires 
a great increase in magnetic field strength even in this case. 

We separately investigated the dependence of the solutions obtained on the plasma- 
occupied fireball radius Rq. In contrast to the analytical model where the plasma tem- 
perature did not depend on this parameter but was determined only by the giant-fiare 
duration, our numerical calculation showed that the temperature of the plasma-occupied 
region decreases with its increasing sizes. The dependence of the magnetic field strength 
on Rq remains similar to the analytical model, but the law 60 ~ Rio is replaced by a 
slightly faster decrease in strength. Figure |5] presents the numerically calculated global 
minimum of themagnetic field strength bmin = ^o™"^'?) at i?o = 5 , 10 and 15 km. It is 
well fitted by the formulas 

SGR 0526 - 66 : bmin ^ 220 R:^q' + 140; SGR 1806 - 20 : bmin ^ 120 i^i^f + 60, 

which are valid both for the isothermal case (7 = 0) and for the case of 7 = —1/2. As 
we see from the plots, the minimum magnetic field strength near i?o — 5 km exceeds 
considerably ij^™") 10^^ G. Therefore, the situation where the plasma occupies a 
fairly extended region Rq > 10 km is of greatest interest in our analysis. In this case, 
the minimum magnetic field strength at the fireball center and the upper bound on the 
field strength obtained from our estimate of the magnetar magnetic dipole losses can be 
represented as 

SGR 0526 - 66 : 5^^"^ ~ 2 x 10^*^ G, Bmd ^ 2 R^^} x 10^^ G; (44) 
SGR 1806 - 20 : eJ™"^ ~ R];^ x 10^^ G, Bmd ^ (2 - 6) ^^o^ x 10^^ G; (45) 
SGR 1900 + 14 : 5^™"^ ~ R^^ x 10^^ G, Bmd - (2 - 3) R:^^ x 10^^ G, (46) 

where eJ™"^ for SGR 1900+14 is the same as that for the SGR 1806-20 giant fiare. To 
estimate the strength Bmd, we used the parameters M^s = 1-4M0 and 9 = n/i and the 
following periods and spindown rates: SGR 0526-66: Pio ^ 0.81, P„io ^ 0.65; SGR 1806- 
20: Pio ~ 0.756, P_io ~ 0.8^8; SGR 1900+14: Pio ~ 0.515, P_io ~ 0.5^1.4, (Mereghetti 
2008). As was discussed above, the modulation of the X-ray intensity for known fiares 
by the magnetar spin period leads to the conclusion that the fireball radius Rq cannot 
differ significantly from the neutron star radius Rns- As follows from the estimates f l44p 
- din]), under the condition Rq ~ Rj^s the upper bound on the magnetic field strength 
of the magnetars under consideration turns out to be lower than the minimally possible 
value required for the suppression of neutrino emission. 

CONCLUSIONS 

In this paper, we estimated the energy losses of a nondegenerate relativistic (T > m) 
electron-positron plasma through neutrino emission based on the magnetar model of an 
SGR giant fiare (Thompson and Duncan 1995, 2001). In the absence of a magnetic 
field, the plasma energy losses through neutrino emission were shown to be too large to 
provide the observed energy release at the LT stage of an SGR giant fiare. It follows 
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from our analysis of the neutrino processes considered that in the case of a strongly 
magnetized plasma eB ^ important for the magnetar model, not only the electron- 
positron annihilation into a neutrino pair ([1]) but also the neutrino synchrotron emission 
process ([5]), which is usually neglected, make a major contribution to the neutrino energy 
losses. The plasma neutrino emissivities in these processes were shown to be significant 
even in the case of fairly strong magnetic fields B > 10^^ G. Thus, when the X-ray 
emission of an SGR giant flare is modeled, the plasma energy losses through neutrino 
emission should be properly taken into account, which was not done by Thompson and 
Duncan (1995, 2001). 

To investigate the main features of the fireball neutrino cooling, we considered a simple 
analytical model. It clearly shows that the photon emission at an arbitrary plasma energy 
cannot exceed some maximum value dependent on the size of the plasma occupied region 
and the magnetic field strength in the plasma. Thus, the minimally possible magnetic 
field strength of a magnetar that provides sufficient suppression of its neutrino emission at 
the long tail stage can be found for the observed energy of a giant flare. We numerically 
modeled the neutrino cooling of the SGR 0526-66, SGR 1806-20, and SGR 1900+14 giant 
flares, including all of the neutrino reactions important for this process. The lower bound 
on the magnetic fleld of these objects corresponding to the energy observed in photons at 
the LT stage of SGR giant flares was shown to disagree with the upper bound from our 
estimate of their magnetic dipole losses. Consequently, the magnetar model of an SGR 
giant flare considered here cannot provide the energetics observed at the long tail stage 
in a fairly wide range of parameters. 

We modeled neutrino cooling under the simplifying assumptions that the flreball tem- 
perature and sizes did not change in the time of an SGR giant flare. Allowance for 
the evolution of these characteristics must lead to a reduction in the energy losses of the 
medium through neutrino emission compared to the model considered. However, it is hard 
to expect that the minimally possible magnetic fleld strength required for the suppression 
of the flreball neutrino emission can be signiflcantly lower than the estimates obtained 
here and can become equal to the upper bounds following from our estimate of the mag- 
netic dipole losses for magnetars. Note also that the magnetar model is not without other 
contradictions either. In particular, an important problem (Malov and Machabeli 2006) 
that arises in attempting to explain the existence of radio emission detected from SGRs 
and AXPs (Malofeev et al. 2005) should be pointed out. 

In conclusion, note that the magnetar model considered here was further developed 
by Lyutikov (2006), Beloborodov and Thompson (2007), and Beloborodov (2009). These 
authors used the more realistic magnetohydrodynamic approach to describe the magnetar 
corona and revealed the effect of additional plasma heating due to the magnetic fleld 
energy. Using this approach to describe SGR giant flares at the LT stage could partially 
solve the problem of energy deflciency discussed here. 
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Figure 1: Scheme of the region of a ball with radius Rq filled with an electron-positron 
plasma and trapped by a poloidal magnetic field with field lines closed on the magnetar 
crust. 
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Figure 2: Dependences of the temperature to(fl) and magnetic field strength ho{h) at the 
fireball center on parameter corresponding to the analytical solution. The lines are drawn 
for Rq = 10 km and the parameters 7 = —1/2 and (3 = —3. 
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Figure 3: Temperature to at the fireball center versus parameter for the SGR 0526- 
66 (a) and SGR 1806-20 (b) flares at Rq = 10 km and (3 = —3. The solid, dashed, and 
dashdotted lines correspond to 7 = —1/2, 7 = 0, and the analytical solution at 7 = —1/2, 
respectively. 




Figure 4: Magnetic fleld strength bO at the flreball center versus parameter for the SGR 
0526-66 (a) and SGR 1806-20 (b) flares at Rq = 10 km and (3 = -3. The solid, dashed, 
and dashdotted lines correspond to 7 = —1/2, 7 = 0, and the analytical solution at 
7 = —1/2, respectively. 
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Figure 5: Minimum magnetic field strength 6™" versus radius -Rio of tlie plasma-occupied 
region for the SGR 0526-66 (a) and SGR 1806-20 (b) flares. The filled circles correspond 
to the calculated values of 6™" at Rq = 5, 10, 15 km and 7 = —1/2, f3 = O.The solid and 
dashed lines represent the fits to this dependence at 7 = —1/2 and 7 = 0, respectively. 
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